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Abstract. We present a model of the growth of water ice on silicate grains in the circumstellar envelopes of 
Asymptotic Giant Branch (AGB) stars and Red Super Giants. We consider the growth of ice by gas grain 
collisions, the thermal evaporation of ice from a grain, and sputtering. Our model contains several improvements 
compared to earlier models, including a detailed treatment of the effects of sputtering, a detailed calculation of the 
radiation pressure on the grain, and the treatment of subsonic drift velocities. In terms of drift velocity between 
the grains and gas in the envelope, we find that the ice formation process can be divided into three regimes: (i) a 
sputtering dominated regime where ice growth is heavily suppressed, (ii) an intermediate regime with moderately 
^ ' efficient condensation and (iii) a thermally dominated regime where ice formation is highly efficient. Sputtering is 

the critical factor which determines if ice formation can occur at all. We find that in Red Supergiants, ice formation 
■ is suppressed, while the winds of OH/IR stars allow for efficient condensation and can convert significant fractions 

of the available water vapor (tens of percent) into ice mantles on grains. Population II AGB stars hardly form 
ice due to their low dust to g as ratios. We also modify an analytical equation describing condensation and 
depletion djura fc Morris! fl985) in order to give reasonable results for high and low drift velocities. Initially, 
ice will condense in crystalline form, but continuing condensation at low temperatures, and damage caused by 
interstellar UV photons favor the production of amorphous ice as well. We predict that a significant fraction of 
' the ice formed will be amorphous. 

Key words, circumstellar matter - infrared: stars - stars: abundances - stars: AGB and post-AGB - stars: mass- 
'_h ' loss 

| 1. Introduction silicates. At the same time, the condensing molecules dis- 
• *h . appear from the gas, and the associated emission of molec- 
^ ■ Asymptotic Giant Branch (AGB) stars represent the fi- ular lines becomes weaker. This applies both to the main 
^ I nal phases of the evolution of stars with low or inter- i ce molecule (e.g. H 2 0) and derived species (e.g. OH maser 
mediate Main Sequence (MS) mass (1 < M < 8 M ). emission). An understanding of the ice condensation pro- 
During the AGB phase, which lasts some 10 6 yr, stars cess is therefore critical for the interpretation of both 
expel matter at a high rate (10~ 7 to 10~ 4 M Q /yr). Mass- the so lid-state and the molecular diagnostics of AGB en- 
loss increases over time, eventually reaching rates in excess velopes 

of 10- 5 M /yr. The outflowing matter creates a dusty Observations reveal that water ice is an important 
molecular circumstellar envelope which may completely solid state compone nt in the circumstellar environments 
obscure the central star. As the dust is driven outwards of ev0 lved stars. Water ice is found in the environments 
by radiation pressure of the central star, it drags the of AGB star s like OH 26.5+0 .6, OH 127.8+0.0, OH 
molecules in the envelope along with it. Given the right 231.8+4.2 (see lOmont et al.ll99fl and references therein) , 
conditions, some species of these molecules may condense post _ A GB stars like HD 161796 faoogzaad et aLlEool 
onto the grains as ice further out in the outflow. The ice planetary nebulae like NGC 6302 (see e.g. iMolster et alJ 
on these grains forms an important diagnostics. It can be Eool. and at least one Red Super Giant, NML Cyg 
detected by its infrared features and provides information jMolster et aiEooi . In these environments it has mainly 
about the physical conditions in the envelope. Ice mantles been observed m a crystalline state. This has to a large ex- 
may also be chemically important, for example for the tend been deduced from the pre sence of both a 43 ^m and 
formation of carbonates ( |Kemper et aLj|2002j) or hydrous 6 Q urn feature in t heir Infrared Space Observatory (ISO, 
iKessler et al .119961) spe ctra. Detailed s tudies of these spec- 
Send offprint requests to: C. Dijkstra (dijkstra@astro.uva.nl) tra can be found in e.g. lBarlowl (|l998) (who describes the 
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featur es for different types of sources) and lSvlvester et all 
l|l999|) (who describe the features for AGB OH/IR stars 
specifically). 

In this paper we present a detailed model which de- 
scribes the growth of water ice on silicate grains in the 
circumstellar envelopes of oxygen-rich AGB stars. The 
model is also applicable to Red Super Giants. The gen- 
eral conditions for the formation of water ice in cir- 
cu mstellar envelopes hav e been described and discussed 
bv lJura fe Morrid l)l985|) : (hereafter refered to as J&M). 
Their model includes most of the major processes involved 
in the formation of ice. However, as was already pointed 
out by J&M, two ingredients are missing in their model: 
sputtering (the mechanical removal of ice lattice parti- 
cles from the surface of the grain by energetic collisions 
with gas particles) and subsonic drift velocities (i.e. sub- 
sonic velocity differences between the dust and gas). We 
will show, that a detailed treatment of these ingredients is 
necessary if we wish to properly describe the ice formation 
process around evolved stars in a broad set of circumstel- 
lar environments. Compared to the J&M model, our model 
also improves the calculation of the radiation pressure on 
the grain, taking into account the optical properties of the 
grain and the spectral energy distribution of the driving 
radiation. 

In Sect. El we give a complete description of the model, 
while in Sect. 01 the results of the model will be presented 
and discussed. In Sect.^Jwe apply our model to an OH/IR 
star, a Red Super Giant and a population II AGB star and 
discuss why ice does or does not grow in the circumstellar 
environments of these stars. The conclusions of this paper 
are given in Sect. 03 Sect. will include discussions on 1) 
the effects of mass-loss rate, grain size and drift velocity 
on ice formation, 2) the distance dependence of ice forma- 
tion and the question whether crystalline or amorphous 
ice will form, and 3) the derivation of a generally valid an- 
alytical expression for the depletion of water vapour from 
the envelope. 

In a future paper (Dijkstra et al., in prep.) we will 
describe the observational applications of our model cal- 
culations in detail. In that paper we will use the model to 
predict the infrared spectra emitted by evolved stars under 
the influence of the ice formation process. The resulting 
spectra will be compared with observations of individual 
stars, and used to derive their stellar parameters. 



2. The model 

2.1. Radiation force, drag force and drift velocity 

A silicate grain near an AGB star, will be blown away from 
the star by radiation pressure. As the grain moves away 
from the star it travels through the circumstellar envelope 
which surrounds the AGB star. On its journey the grain 
will hit gas molecules in the envelope. The collisions create 
a drag force on the grain. We will assume that the drag 
force will be in equilibrium with the radiation force at any 



time and place. We define r to be the distance of the grain 
from the center of the star. 

The radiation force on the grain is given by 



-^rad 



47rcr 



/>oo 

' / Qabs.A-^AdA 
J A=0 



(1) 



where Q a bs.A is the absorptivity efficiency of the grain 
at wavelength A, and H^ is the local Eddington flux at 
wavelength A. For the calculation of H\ we use th e radia- 
tive transfer programme MODUST (|Bouwmanl200lf) . which 
models the (emergent) spectrum of a spherical dust shell 
around a central star, a is the radius of the grain and c is 
the speed of light. Eq. ^may be simplified to 



-^rad 



4c 



(2) 



where L* is the total luminosity of the star, i.e. the lumi- 
nosity integrated over all A. Qabs is the mean absorptivity 
efficiency integrated over all A with respect to the local 
energy distribution. In Sect. l2~7l the calculation of Q a bs,A 
will be discussed in detail. 
The drag force is given by 



fdrag = 7ra 2 /w t 2 h 



64 

97T* 



(3) 



as discussed bv lDominikl lll992). Here p is the mass density 
of the gas in the envelope, while i>th is its thermal speed, s 
is defined to be v^m/vth where thrift is the drift velocity, 
i.e. the velocity of the dust grain with respect to the gas. 

Eq. ^ or 121 together with Eq. can be solved for the 
drift velocity of the grain. Setting F ra( j equal to i^rag we 
obtain 



"drift = Vth 



\ 




(4) 



"drift is an important quantity for the model since, as 
we will show later, for supersonic drift velocities it de- 
termines both the rate at which water molecules are col- 
lected from the gas and the rate at which these molecules 
are sputtered from the grain's surface again, thus influ- 
encing the amount of ice that forms. Wdrift also enters the 
expression for da/dr, the rate of change in radius of the 
icemantle as a function of distance to the star, which we 
will derive now. 

If dV is the change in volume of the grain as it collects 
dNji 2 o water molecules, each occupying a volume Vh 2 o, 
from the gas in the envelope, then we have 



da 



dV dN H2 oV H2 o 



Ana 1 



Ana 2 



(5) 



The above equation assumes that all water molecules are 
distributed uniformly over the grain surface, i.e. the grain 
remains spherical at all times. For a given time interval dt, 
dNn 2 o is given by (i) the rate at which water molecules 
from the gas in the envelope are collected by the grain, i.e. 



C. Dijkstra et al.: Water ice growth around evolved stars 



3 



the growth rate Z gT , (ii) the rate at which water molecules 
are lost from the grain due to evaporation, i.e. the evapo- 
ration rate Z CV1 and (iii) the rate at which water molecules 
are lost from the grain due to sputtering, i.e. the sputter- 
ing rate Z sp . With the velocity of the dust particle with 
respect to the star given by v^ust = dr/dt and the density 
of the ice in the mantle given by p lcc — y/"? , where m p 
is the mass of a proton, we find 



da = 



(Z gT — Z cv — Z S p)dtVn 2 



Ana 2 

(Z gI — Z cv — Z sp )Vn 2 c 



Ana 2 



■dr, 



^dust 



(6) 



yielding 



will happen at low temperatures and with water molecules 
sticking to an ice surface. Because of the dipole forces in 
water ice, sticking is likely in this case. However, initially 
ice will grow on silicate cores and possibly on other ma- 
terials, where sticking has to rely on van der Waals forces 
only. In this case it is possible that a will drop below 
unity. Also at high drift velocites sticking may become in- 
efficient. While we don't treat the velocity dependence of 
a, ice growth becomes automatically inefficient at those 
drift speeds because of sputtering. 

For the limiting cases of small and large drift velocities 
Z gr may be written as 



I 9 V 2 

Z T =\ a4?ra n H 2 o^^ : s H . 2 o < 1 
I cnra 2 nH,o|vdrift| : «h 2 o > 1 



(12) 



9m,, 



da 



dr 2np- lce v Aust d- 



Zr. 



(7) 



v drift enters Eq. [7] if we note that it can be expressed 
as Adrift = «dust — Vgas, where i> gas is the speed of the gas in 
the envelope with respect to the star. v gas is assumed to 
be an input parameter of the model. With this expression 

and Z sp in more detail. 



for ^ we will now study Z gr , Z c 



2.2. The growth rate, Z gT 

Assuming a Maxwell velocity distribution for the gas, the 
rate at which the d ust grain collect s water molecules from 
the gas is given by jDomini]J ll992) 



_ ana 2 v^° n H2 p 

Z gr — — X 

V 71 " S H 2 Q 



^(2s| 2O + l)0( SH2 o) + SH 2 oe 



(8) 



where iih 2 o is the number density of water molecules in 
the gas, a is the sticking probability of a water molecule 
on the grain's surface. 



sh 2 o 



Adrift 



(s H2 o) = —j= \ e~* 
V"" Jo 



dt 



(9) 
(10) 



2.3. The evaporation rate, Z cv 

Consider a mass of water ice surrounded by water vapour. 
The vapour pressure of water, p^ap 1 is defined to be the 
equilibrium pressure at which the rate in which the water 
ice evaporates from a surface is equal to the rate in which 
the water vapour condenses into ice on the surface. If the 
water gas pressure is larger (lower) than p^ap ice will grow 
(evaporate). We can (based on the Clausius-Clapeyron re- 
lation) express p^ap as 



Pvap° =PH 2 



vap 

,H 2 
''vap 



™H 2 



1.6 x 10 12 e 



(13) 



where ph 2 o is the water gas pressure in dyn cm 



,h 2 o 

vap 



IS 



the number density of water vapour molecules at phoO 
it. is the 

aadll2()()ll). We assume that Td us t = T eas . We 



PHoO = 

l|Kouchil 



P^rP and TH„st,_is the temperature of the grain 
198^ IS" 

will need p^* p ° in our calculation of the evaporation rate, 

Z C v 

The evaporation rate can be calculated from p^£ p > us- 
ing that in thermal equilibrium (TE) the growth rate must 
equal the evaporation rate, i.e . Z„ = Z „ Y . The growth rate 



for a grain in TE is given bv iDominikl l)l992() as (see also 
Ea.[T31in the limit s H2 o < 1) 



a4™ 2 n« 2 p 



h 2 o 
2y^ 



Using Z e 
we find 



Z„ 



a 



(14) 

1 and combining Eas. 1131 and 1141 



is the error function and 



h 2 o 



2fcTga 

18m r 



(11) 



is the thermal speed of the water molecules in the gas. 
Tga S , the temperature structure of the gas and dust in 
the envelope as a function of distance to the central star, 
may be given as a power law or calculated using radiative 
transfer codes. 

The sticking probability will be assumed in our calcu- 
lations to be equal to one, since most of the ice formation 



47ra 



2 h 2 o 

vap 



W H 2 

20F 

^HaO H 2 



47TO 71h 2 



vap 



th 



Ph 2 o 20? 



H 2 



= 2™ 



2 v th_ n H 2 Q 
7T S 



(15) 



where we have defined S = Ph 2 o/p^ p as the supersatu- 
ration ratio. If S < 1 ice evaporates from the surface, if 
S > 1 it will grow (as discussed at the beginning of this 
section). 
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2.4. The sputtering rate, Z sp 

As the grain moves outwards it experiences head-on col- 
lisions with gas particles in the envelope. Most of the gas 
in the envelope is present in molecular form and its main 
constituents are molecular Hydrogen (H 2 ), Helium (He), 
Carbon monoxide (CO) and molecular Nitrogen (N 2 ). As 
these particles collide with the grain, they may, for large 
enough drift velocities, mechanically remove lattice par- 
ticles from the ice mantle of the grain. This is known 
as sputtering. We assume that the original silicate grain, 
which is the core of the dust particle, is not destroyed by 
sputtering 1 . To incorporate the effec ts of sputtering in ou r 
model, we follow the treatment bv IWoitke et alJ lfl9 93b 
who describe the sputtering process in detail. We have 
used their equations and will write them down below for 
completeness. We ignore chemisputtering, i.e. the removal 
of ice lattice particles by means of chemical reactions with 
other species. 

In order to calculate the sputtering rate for a species 
i, Z sPt i, we first derive the collision rate of the dust grain 
with the species. In the regime of supersonic drift velocities 
this rate is given by 



ry 2 

^coi,i = riiiia Udrift 



(16) 



where rii is the gas number density of species i. The as- 
sumption of supersonic drift velocities at this point does 
not affect the results of our model, since in the ice forming 
regions (as we will show later) sputtering only occurs at 
supersonic drift velocities. Second, we calculate the sput- 
tering rate for each species using 



Zsp,i — Z C xA.il£ 



co\.i * sp,? 



where Y sPt i is the sputtering yield f or species i. Y sri _j can be 
calculated from llRoth et al.lll978irBohdanskv et al.lll98(A 
IStrazzulla et al.lll985[) " 



Y„ 



0.0064** 7i s /3( 







_ £c,-,, N2 
Ei > 

E% ^> E cr i 

E < E 

s - J % _ J -'cr,« 



(18) 



In these equations Ei = 0.5miVdrift 2 is the impact ki- 
netic energy of the species i as its collides with the grain 
7 . mi is the mass of an individual mem- 



and 7^ — 

ber of a species i, and mo is the mass of a (target) lattice 
particle. In our case this is equal to the mass of a wa- 
ter molecule, i.e. mo = 18m p . The sputtering threshold 
energy, £ crii , is given by 



I 



E — ) 7<(l-7*) 



1 8£ b (^r 



-<0.3 

J 

- > 0.3 



(19) 



1 The silicate core is held together by chemical bonds, while 
the ice mantle is held together by physical bonds. Chemical 
bonds are generally stronger than physical bonds. This makes 
them less sensitive to sputtering. 



where Ey, is the surface binding energy of a lattice par- 
ticle on the target material. Sputtering shows a thresh- 
old behaviour. Only particles with sufficient energy, i.e. 



Ej, > E r 



are able to remove lattice particles. Particles 



with E. t < E cl i are not. The value of E^ differs for each 
material and is lower for water ice (£h = 0.1 ... 0.37 eV) 
than for silicates like MgSi0 3 (E h = 4.5 eV). Notice that 
£b is roughly an order of magnitude larger for MgSi03 
than water ice, which at least qualitatively justifies our 
assumption above that only the icemantle is affected by 
sputtering and not the silicate core. Representative values 
of i?b for different speci es, including those m entioned here, 
are listed in Table 3 of IWoitke etafl l|l993|) . 

Finally, with the sputtering yield for each species cal- 
culated we can now calculate the sputtering rate for each 
individual species and from this the total sputtering rate 
Z sp using (if there are N species) 



Z S p 



N 



(20) 



Another process which may lead to the removal of ice 
from the grains is grain-grain collisions. Grains of differ- 
ent sizes and/or optical properties will drift through the 
enevelope at different speeds and can therefore collide at 
significant speeds. Since the current models assume a sin- 
gle grain size, relative velocities between grains are zero 
and this effect can be ignored. However, in a future study 
with grain size distributions, it should be taken into ac- 
count. 



(17) 2.5. Depletion 



As the grain moves outwards it collects water molecules 
from the gas, which are then no longer available to interact 
with new grains, i.e. the water molecules are depleted from 
the gas. The effects of depletion are incorporated into our 
models, by calculating the effects of the passage of grains 
on the number density of water molecules in the gas, nu 2 o- 
We have 



dr 



dt 



2"h 2 o 



(21) 



for the change in nu 2 o as grains move a distance dr away 
from the star in a time interval dt and collect A/Vh 2 o water 
molecules from the gas. ridust is the number density of dust 
grains in the envelope. The first term on the right hand 
side of Eg. 12 II takes into account the depletion effect. The 
second term takes into account the effect that as grains 
move outward, the surrounding gas dillutes into a larger 
volume and therefore iih 2 o drops. Dissociation of H 2 
molecules by the interstellar UV field may also lead to a 
decrease of the available molecules. We have ignored this 
effect. It will only be important for low mass loss rate 
stars, where (as we will show below) sputtering prohibits 
ice formation. 
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"■dust is calculated using 



"■dust 



3M/ dg 



167r 2 a 3 p si ir 2 i; dust 



(22) 



where M is the gas mass-loss rate, /d g is the dust-to-gas 
ratio, ao is the radius of the silicate grain, and hence the 
initial radius of the core mantle grain that forms as ice 
grows on it, and p s n is the bulk density of a silicate grain 
(which we took to be 3300 kg m~ 3 ). Due to the generally 
small values for /d g (/dg~0.01) we assume that M is not 
just the gas mass-loss rate, but the total mass-loss rate. 

2.6. Crystalline or amorphous ice 

Ice can exist in a crystalline or amorphous state. 
Crystalline ice can form if the mobility of the water 
molecules on the grain surface is large, and the flux of 
water molecules on the surface is relatively low. In this 
case, each water molecule has sufficient time to relax to a 
low energy configuration and hence help in the creation of 
a crystalline lattice. Amorphous ice forms when the mo- 
bility of the water molecules on the grain surface is small 
and the flux of water molecules on the surface is relatively 
high. Now the water molecules have no time to relax to a 
low energy configuration, resulting in an amorphous struc- 
ture. A high (low) mobility of the water molecules can be 
achieved at high (low) temperatures. 

In the infrared, crystalline ice will emit soli d state 
features at 43 /jm and 60 fim l|Bertie et al.lll969l) . while 
amorphous ic e will emit at 43 /^m but not at 60 /im 
l|Moorel ir999). In order to correctly predict and/or ex- 
plain the infrared spectra emitted by the circumstellar 
envelopes of evolved stars, and hence the physical and 
chemical conditions inside these envelopes, the question 
whether crystalline or amorphous ice forms must there- 
fore be addressed. We will present the spectra emitted 
by (the grains in) the circumstellar envelopes in a future 
paper. For completeness however, we already include the 
crystalline/amorphous ice discussion in this paper. 

To investigate whether crystalline or amorphous ice 
will form we check if th e following condition, originally 
derived bv lKouchi etaD (|l994|) - is full-filled 



j ^ *^ crit 



(23) 



where T gT is the growth flux of water molecules per unit 
area and JF cr it is the socalled critical flux. T gT is given by 



while Tc-cit is given by 
£) s (0)e" ( ^&r } 



J~ crit — 



^lattice^ 



(24) 



(25) 



where D s (0) = 1.74 x 10 1 m 2 s 1 , k is Boltzmann's con- 
stant and E s /k — 4590 K. <zi a ttice is the lattice constant of 



water ice, which is 4.5 A. We will once more assume that 
Td U st — T gas (see Sect. 12. 3|) . Condition l2~3l summarizes the 
above considerations about whether crystalline or amor- 
phous ice will form. Indeed, if, at a certain temperature 
(and therefore a certain mobility of the water molecules 
on the surface of the grain), the flux of water molecules 
on the surface is lower (larger) than a certain critical flux, 
crystalline (amorphous) ice will form. For completeness, 
we note that a discussion on the formation of crystalline 
or amorphous Carbon, in essence similar to t hat g iven here 
for water ice, is given bv lOail fc Sedlmavrl lll98^ . 

A third mechanism that determines the type of ice that 
is formed, besides fluxes and temperature, is UV radiation 
from the Inter s tellar Med ium (ISM). Experiments d one by 
iLenault et ail l|l983h and lKouchi fc Kurodal l|l990h reveal 
that UV radiation can alter crystalline ice into amorphous 
ice when the temperature is below 70 K. This happens un- 
der the release of H2 gas. Above 70 K the crystalline water 
ice is able to repair itself from UV ra d iation damage and 
remain crystalline. iMoore fc Hudson! l|l992F) have shown 
that the unaltered fraction of crystalline water ice under 
the influence of UV radiation, can be given by 



$(£>) = e 



(26) 



In this equation Q(D) is the unaltered fraction of 
crystalline water ice after it has been irradiated with a 
dose D of UV radiation. D is given in eV. k is given in 
molecules/eV and is a measure for the amount of water 
molecules in the crystalline lattice that can be reordered 
into an amorphous structure for each eV of UV radiation 
supplied, k depends on the temperature of the ice. A graph 
showin g this dependence is given by IMoore fc Hudsorl 
(1992): k increases with decreasing temperature, mean- 
ing that for lower temperatures one eV of UV radiation 
can be used to amorphousize a larger fraction of ice. 

At present, our model does not include the effects of 
the UV radiation from the ISM, since we expect its effects 
on the ice formation process to be negligible. However, we 
do recognize its importance regarding the type of ice that 
is created, and hence the spectrum that the final grain 
will emit. In Sect. [3] we will therefore include the issue of 
the UV radiation from the ISM into our discussion where 
needed. 



2.7. Calculation ofQ^^x 

As the grain collects ice it increases in radius and its opti- 
cal constants, and hence Q a bs,A, will change. The change in 
Qabs.A will affect -F ra d, and so thrift- This in turn will influ- 
ence the (final) radius of the grain. We therefore calculate 
Qabs.A at each r to incorporate this effect. Based on the 
value of the size parameter, x = we either ap p ly the 
routine bhcoat (taken from IBohren fc HuffmannI l|l983h 
: used when x < 30) or an effective medium calculation in 
combi nation with a MIE calculation ijBohren fc HuffmannI 
l(l983h : used when x>30) to determine <3abs,A- The optical 
constants used in our calculations for the silicate core and 
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the ice m antle were taken fro m lSuhl <|1999) (cold silicate 
data) and IBertie et al.l l|l969|) respectively. 

3. Results and discussion 

In this section we present and discuss the results of the ice 
growth model described in the previous section. First, we 
will discuss the effect of M and the initial grain size, ao, on 
the formation of ice in the circumstellar envelope. Second, 
we will describe the effect of i>drift and a on the deple- 
tion of water vapour from the envelope. Third, we discuss 
the distance dependence of ice formation and address the 
question whether crystalline or amorphous ice will form, 
followed by a discussion on the type of spectra produced 
by the core mantle grains formed in the envelope. Fourth 
we derive a generally valid analytical expression for the 
depletion of water vapour from the envelope. Finally, we 
apply our model to individual stars and compare our re- 
sults with observations on these stars. 

Our discussion on the above topics will largely be based 
on models of circumstellar envelopes with M ranging be- 
tween 10 -7 and 10 -3 M Q /yr. The choice of modelling M 
above the more conventional upper limit of 10 -4 M Q /yr is 
based on observations of the post Red Sup er Giant binary 
system AFGL 4106 dMolster et all Il999h and the post- 
AGB stars HD 161796 llHoogzaad et a.l.ll20ol and IRAS 
16342-3814 (Dijkstra et al., in prep.). All these systems 
are thought to have experienced mass-loss rates of sev- 
eral times 10 -4 M Q /yr up to 10" 3 M Q /yr at some point 
during their history, or at least have had densities in 
the wind that, assuming spherically symmetric mass-loss, 
would require mass-loss-rates of the order of 10~ 3 M Q /yr. 
Calculations of models with M above 10 -4 M Q /yr there- 
fore seem justified. The models contain a number of as- 
sumptions. Each envelope is composed of single sized par- 
ticles with radii of cither 0.01, 0.10, 1.00 or 10.0 fim. The 
initial abundance of water vapour (by number with re- 
spect to H) is £q 20 = 1 x 10~ 4 . The dust-to-gas ratio is 
equal to /dg = 0.01. The velocity of the gas is given by 
u gas = 15 km/s. The central star has an effective tempera- 
ture of T e g = 3500 K and a radius R± — 300 R Q , yielding 
a luminosity L*»T.2xl0 4 L Q . We will address the impor- 
tance of these assumptions where needed. 

We used the radiative transfer programme modust 
(see Sect. I2.1fl to provide a correct physical description 
of the temperature and radiative flux in the envelope. 
For (very) high mass-loss rates, and hence high optical 
depths, the modust models do not fully converge. When 
this problem occured, we resolved it by moving the inner 
radius of the circumstellar envelope outwards. Our results 
will not be severely affected by this. 

3.1. Mass-loss rate and initial grain size 

In Fig. n the relative increase in the mass of the core man- 
tle grain (when the complete ice mantle has formed), ^p, 
is shown as a function of M for ao = 0.01, 0.1, 1.0, and 
10.0 /xm. Fig. ^ shows that for a given ao ice formation 



Fig. 1. The relative increase in the mass of the core mantle 
grain (when the complete ice mantle has formed) as a 
function of M for initial grain sizes ao = 0.01, 0.1, 1.0 and 
10.0 /zm. For more details see Sect. 13. ll 
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becomes more efficient as M increases. For an increasing 
M, the density, p, increases and the drift velocity, thrift 
decreases 2 . The increasing density increases the growth 
rate of the ice on the grain, while the decreasing drift ve- 
locity lowers the effect of sputtering on the grains, leaving 
newly formed ice layers intact and promoting ice forma- 
tion as well. Fig. also shows that for smaller values of 
ao, significant ice formation starts at lower M. For smaller 
ad ^drift is lower and sputtering is less important, again 
leaving newly formed ice layers intact. In Sect. 13.21 we will 
come back to the importance of sputtering. 

We conclude that small grains have the best ability 
to grow ice and will most likely dominate the ice forma- 
tion process if present in sufficient relative numbers. In 
grain size distributions w e generally expect sm all parti- 
cles to dominate (see e.g. iDominik et aUll989|) . We may 
thus expect that models with small single sized particles 
and grain size distributions will yield similar results. If we 
for the moment take the behaviour of the smallest parti- 
cles as typical, the curve for ao = 0.01 /im suggests that 
modest ice growth may already start at mass-loss rates 
of 10 -6 Mo/yr and that major ice formation only sets in 
when the mass-loss rate is typically of the order of several 
times 10 -5 to 10 -4 M©/yr, i.e. during the SW phase of 
an AGB star (see Sect.^l. 

3.2. Mass-loss rate, drift velocity, initial grain size and 
abundance 

Fig. |2 shows two contour plots of /(oo), the final water 
vapour abundance in the circumstellar envelope after ice 



2 As M and hence p increaeses the drag force on the grain 
increases. The drag force must balance the radiative force on 
the grain, which remains approximately constant. Therefore 
«drift must decrease. 
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Fig. 2. Contours of constant final water vapour abundance relative to the initial abundance, /(oo), as a function of ao 
and M (upper panel) and ao and Vdrift (lower panel), /(oo) is given in percentages. Also shown are lines that indicate 
the critical velocities for sputtering by H2, He, CO and N2, and the sound velocity at 100 K. 
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formation relative to the initial abundance, £q • In the 
upper panel of the figure /(oo) is shown as a function of 
do and M and given in percentages. The contours, except 
for the 99% contour, are shown at 5% intervals. The unre- 
alistic small scale curvatures in some of the contours (see 
e.g. the 90% contour) are due to the interpolation of the 
contour values on the low resolution grid and may be ig- 
nored. M is varied between 10 -7 M Q /yr and 10 -3 M Q /yr. 
These values correspond to those in Fig.^ The lower panel 
shows the same contour plot, but this time with Adrift in- 
stead of M on the horizontal axis. thrift decreases as M 
increases. Also shown in both contour plots are lines indi- 
cating where Adrift reaches the critical velocity for sputter- 
ing by H 2 , w critH2 (~9 km/s), He, -u critHc (-6 km/s), and 
CO and N 2 , w critco = v clitN2 (~3 km/s), and the sound 
velocity at 100 K, c s (~0.6 km/s). The 'No Ice Formation' 
labels in these plots identify the regions where no ice for- 
mation occurs. The 'No Stars' label in the lower contour 
plot indicates the region where M exceeds 10 -3 M /yr 
and no stars are expected observationally. 

The upper panel of Fig. shows that /(oo) decreases 
as M increases and/or ao decreases. Also, ice formation 
becomes more and more efficient, as is suggested by the 
increasing density of the contours. As was discussed in 
Sect. 13.21 an increase in M leads to higher densities and 
lower drift velocities, and hence larger growth rates and 
less sputtering. This promotes the formation of ice and 
makes /(oo) smaller. A decrease in ao also leads to lower 
drift velocities, i.e. less sputtering, and (at a fixed M) it in- 
creases the total collecting area of the grains. Again, this 
promotes the formation of ice and makes /(oo) smaller. 
The obtained values for /(oo) cover the full range be- 
tween 100 and 0%. /(oo) = 0%, which corresponds to 
a total freeze out of all water vapour from the circumstel- 
lar envelope, is obtained in the most extreme cases, when 
M«10" 3 M /yr and a «0.01 fim. 

The contours of /(oo) in the upper panel of Fig.|2run 
parallel to the lines which indicate where the drift velocity 
reaches the critical velocity for sputtering by H2, He, and 
CO and N2, suggesting that the drift velocity is (indeed) 
a very important quantity in determining how much ice 
will grow. This is illustrated in the lower panel of Fig. 

In the lower panel the contours depend (almost) only 
on Adrift and not on ao- This allows us to characterize the 
ice formation process in terms of Vdrift ■ The lower panel of 
Fig-Elshows that for any given particle size, ice formation 
is impossible when u drift > w C rit H2 ■ When t> dr ift < v critll2 
ice formation is switched on, even though He, CO and 
N2 still cause sputtering. The amount of ice formed in 
this regime is modest however (/(oo) > 92%). Sputtering 
ceases below u d rift = ^ciitco = u crit N ■ From now on we 
will define the regime above w dr ift = iW co = u crit N2 as 
the sputtering dominated regime. Between Wdrift~c s and 
Adrift = v CT it co = v C rit N2 there is an intermediate regime 
where ice formation becomes increasingly more efficient 
(as can be seen by the increasing contour density), but 
still does not exceed 20%. At subsonic drift velocities, 
i.e. when thrift < c s , /(oo) seems to become very sen- 



sitive to Adrift, and a total freeze out of water vapour 
from the envelope occurs for the lowest drift velocities 
(vdrift~6xl0~ 3 km/s). 

We emphasize that in the subsonic or thermally dom- 
inated regime, it actually is not u d rift that controls the 
ice growth, but it is M that does. Eq. ^| shows that for 
subsonic drift velocities the growth rate of ice on a grain 
does depend on the density p, and thus M, but not on 
Adrift- Only the thermal motions and density of the gas 
are important in this regime. The apparent rapid decline 
of /(oo) with thrift i s thus simply the effect of the increas- 
ing mass-loss rate of the star. In the sputtering dominated 
regime it is also M that we change, but now it is indeed 
the effect of v drift that is truely important, since here it 
controls the sputtering yields (see Eq. IT%|l , and thus for 
which gas species sputtering is switched on or off. In the 
intermediate regime, both M and thrift are important. In 
this regime M sets the density in the envelope, and thus 
the growth rate, and it partially sets the value for i>drift 
(^drift depends on more variables than M alone) . Udrift also 
sets the growth rate, together with the time that the dust 
spends in a certain (ice formation) region of the envelope. 

Fig. HI can be used to derive the abundance of water 
vapour in the envelope of an evolved star, keeping in mind 
the assumptions discussed in Sect. |3 The assumptions 
that the dust-to-gas ratio is given by / dg = 0.01 and the 
initial abundance of water vapour by e^ 20 = 1 x 10~ 4 are 
an important aspect of Fig. [2J When /d g is lower (larger) 
than /d g = 0.01, the figure only provides an lower (upper) 
limit on /(oo), since there will be less (more) dust particles 
available to grow ice on. The effect of a changing e^ 2 ° is 
not trivial since it is dependent on multiple variables that 
continuously change throughout the envelope (see Ea. l21|) . 
Our models suggest that when ep 12 is lower (larger) than 
£q 2 ° = 1 x 10~ 4 , Fig. [3 provides an lower (upper) limit 
on /(oo). This is due to the overall decrease (increase) in 
the partial pressure of the water vapour, which brings the 
starting point of ice formation more outwards (inwards) 
in the envelope, and the decrease (increase) of the total 
collecting area of the grains at any given distance. In Sect. 
13.41 an analytical equation for /(oo) will be derived that 
is generally applicable. 

3.3. Crystalline and/or amorphous ice? 

We have investigated the question whether crystalline or 
amorphous ice will form in the circumstellar envelope. 
Our model shows that both types of ice may be present. 
As a typical example we briefly consider a star with 
M = 10~ 5 M Q /yr and a = 0.01 fxm. We find that as the 
grain moves away from the central star, initially crystalline 
ice will be formed, meaning that condition 1231 is met. At 
a distance of approximately 5100 AU and onwards, amor- 
phous ice is formed. Since in condition 1231 JF gr decreases 
as a function of distance to the star, it must be the de- 
crease in Tdust which causes the breakdown of this condi- 
tion. In this example, Td us t ~ 65 K when the transition 
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occurs. So from 5100 AU and onwards, each new water 
molecule that lands on the surface of the grain is immedi- 
ately immobilised, preventing crystalline ice to form. We 
find that (not shown) the grain has reached about 70% 
of its final radius, at the time the transition occurs. This 
means that a substantial fraction of the grain will be com- 
posed of amorphous ice, in this case a volume fraction of 
100(1 - (70/100) 3 ) w 66%. 

Fig. |31 shows the dependence of ice formation on the 
distance towards the central star. For a given initial grain 
size, the figure shows, as a function of mass-loss rate: 

• -Rstart, the distance where ice formation starts (star 
symbols) 

• -Rstop, the distance where ice condensation terminates 
because of too low H2O densities (diamond symbols). 

■ -R70K, where the temperature drops below 70K. At 
higher temperatures, UV damage in crystalline ice is 
immediately repaired, while at lower temperatures, UV 
radiation converts crystalline ice into amourphous ice 
(squares) . 

. i?CrAm, the distance outside which ice already con- 
denses in amorphous form (triangles). 

. i?uV: the distance to which interstellar UV radiation 
can penetrate (crosses). Outside this radius, ice will be 
amorphisized unless the temperature is above 70K. 

We have defined i? s t op to be the distance where the grain 
has reached 98% of its final size. We determine the dis- 
tance to which the UV radiation can penetrate into the 
envelope, Rjjy, by calculating where the UV optical depth 
of the envelope, t\jv (at wavelengths A < 3900A), reaches 
unity. Generally, at a distance r from the central star we 
have (using Eq. 



T\jv(r)= 
and thus 



Qvv^a n dust {r )dr = z (27) 



167ra p s iiVdustr 



3Quv/d g Af 



167ra /3 s iiWdustTuv(^) 



(28) 



From this we calculate i?uv= r ( T uv = !)• We assume that 
Quv, the absorptivity efficiency for UV radiation by the 
grain, is given by Quv = 1 and v gas = 15 km/s. The 
shaded areas in Fig. [3] indicate those regions where crys- 
talline ice is formed. In the unshaded regions amorphous 
ice forms. 

Fig. shows that increasing the mass-loss rate gener- 
ally moves iistart, i?sto P , -R70K and RcrAm inwards. So, ice 
formation starts closer to the central star, the final particle 
size will be reached earlier, and the transition point where 
amorphous ice instead of crystalline ice will be formed 
moves inwards. Meanwhile, i?uv moves outwards, due to 
the increasing optical depth of the circumstellar envelope. 
For envelopes with low optical depths (i.e. envelopes with 
low mass- loss rates and or large dust particles), the ef- 
fect of the UV radiation is generally important through- 
out most of the envelope, and will most likely contribute 



to the formation of amorphous ice. However, the effect of 
the UV radiation can be reversed in regions close to the 
central star, where the temperature exceeds 70 K. 

An increase of the initial grain size generally has a 
modest effect on the distance dependence of ice formation. 
Comparing the two panels in Fig. |3 it can be seen that 
-Rstart, -Rstop, -R70K and -R&Am only show minor changes 
in general. However, i?uv changes over an order of mag- 
nitude, getting closer to the star as ao increases. This is 
again due to the decrease in the optical depth of the enve- 
lope. The effect of UV radiation from the ISM on the type 
of ice that forms is thus most important for stars with on 
average large particles in their envelopes. 

As an AGB star evolves its mass-loss rate increases 
over time. Fig. then suggests that the nature of the ice 
formation process will change over time as well. So, to get 
a qualitative impression on the evolution of the distance 
dependence of ice formation, one might wish to substitute 
the mass-loss rate on the horizonatal axis by a measure 
for the time that the AGB stars has spent on the AGB. 
It must be noted however, that as the star evolves, its 
luminosity will increase as well, an effect that is not in- 
corporated in Fig. 

To illustrate the use of Fig. |3 we will look at an ex- 
ample given by the dotted line in the left panel. This ex- 
ample shows how the distance dependence of the ice for- 
mation process can be studied for a star which loses mass 
at a rate of M = 10 -5 M Q /yr and has particles of ra- 
dius ao = 0.01 /im in its envelope. As a grain moves away 
from this star, the first ice will condense on its surface at 
~1500 AU. This ice will be crystalline. In this model we 
found that the temperature at the start of ice formation 
is approximately 94 K (not shown). For the other models 
we typically we found values between 92 and 104 K. At 
^5000 AU the temperature becomes too low to form crys- 
talline ice. Instead amorphous ice will form. Finally, from 
^20000 AU onwards, interstellar UV photons will (start 
to) bombard the surface of the grain, most likely destroy- 
ing all of the crystalline ice that is present. At ^40000 AU 
most of the ice will have formed on the surface. Notice 
that if we had taken a star with M = 10~ 4 M /yr, the 
UV photons of the ISM would only have become impor- 
tant after ice formation had already ceased. Again, any 
crystalline ice still present in the grain would most likely 
be destroyed. 

3.4. A recipe for making ice 

J&M derived an equation for the final abundance, / (00) 
(also see Sect. I3.2JI . of a condensable species, like water, 
in the outflow of an evolved star. The model of J&M con- 
siders the outer circumstellar envelope of the star, where 
the matererial has reached its terminal outflow velocity, 
i> gas . It is assumed that there is a condensation radius, ro, 
such that for r < ro, there is no condensation, while for 
r > ro almost every molecule sticks to the grain. Ice is 
collected by the grains through gas grain collisions and it 
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Fig. 3. Crystalline and amorphous ice formation, for stellar envelopes with 0.01/im(left panel) or 0.1/im(right panel) 
grains. Wind trajectories run as vertical lines through this diagram from bottom to top. The different lines show 
locations where conditions for ice formation change (see text). The shaded area indicates where crystalline ice is 
formed and can survive. 




is removed by thermal evaporation. The depletion of wa- 
ter from the gas phase is taken into account as well. Two 
very important assumptions of the J&M model are that 
1) sputtering is unimportant, since J&M expect the drift 
velocities of the grains to be less than 5 km/s, and that 2) 
the grains are moving through the gas at supersonic drift 
velocities. Based on our results in Sect. 13.21 we can predict 
that assumptions 1) and 2) will cause the J&M model 
to break down in respectively the sputtering dominated 
and thermally dominated regime. In the sputtering domi- 
nated regime, J&M will predict ice growth, while sputter- 
ing will heavily supresses or completely block ice growth. 
In the thermally dominated regime the growth rate of the 
ice (also see Eq. Il2")) will be underestimated, and hence 
the amount of ice that grows. The equation of J&M will 
only be valid in the intermediate regime, where assump- 
tions 1) and 2) do apply. In this section we will study the 
equation given by J&M, by applying it to the circumstel- 



lar envelopes described at the beginning of Sect. We 
will compare its predicted depletion for these envelopes to 
the depletion predicted by our model calculations of these 
envelopes. Based on this comparison, we will provide a 
modified version of the J&M equation that is valid in all 
regimes, i.e. the sputtering dominated, intermediate and 
thermally dominated regime. 

According to J&M the fraction of a species that re- 
mains in the gas phase, /(oo), is given by 

/(oo) = exp(-a(« dri ft/ug as )Tuv(n))/Quv) (29) 

where Qtjv ~ 1 and ttjv(^o) is a measure of the UV optical 
depth between ro and oo. All other symbols have their 
usual meaning. The depletion of the species from the gas 
phase is (in general) given by 

A = 1 - /(oo) (30) 
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Fig. 4. (Left panel) The depletion (in percentages), as originally predicted by J&M, Ajm, against the depletion 
calculated by our model, A. The solid line shows where Ajm = A. The inset shows a close up of the plot near the 
origin. (Right panel) Same as the left panel, but this time the constraints from Eq. 03 are forced upon Eq. |2H1 (or 
better, Ea. 13 For details and a dicussion see Sect. 13.41 
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In the left panel of Fig. 0] we plot the depletion as pre- 
dicted by J&M, Ajm, against the depletion predicted by 
our model, A, for the circumstellar envelopes described at 
the beginning of Sect. Ajm (plotted on the horizontal 
axis) and A (plotted on the vertical axis) are both given 
in percentages. The solid line shows where A = Ajm- The 
inset shows a close up view of the plot near the origin. It 
can be seen that for A > 20% the J&M model underes- 
timates the depletion, i.e. it overestimates the true /(oo). 
The points for which this happens are indicated with stars. 
Also, for A = 0, the J&M model predicts that the deple- 
tion is larger than zero, i.e. Ajm > 0. Here J&M thus 
overestimate A and underestimate the true /(oo). The 
points for which this happens are indicated with crosses. 
Finally, for < A < 20%, the predictions of the J&M 
model generally agree with our calculations (these points 
are indicated with diamonds). For small/large values of A 
in this range, J&M predict somewhat larger/lower values 
for the depletion however. 

For reasons discussed at the begining of this sec- 
tion, the disagreement between the J&M model and our 
model indeed occurs in the sputtering dominated regime 
(crosses) and the thermally dominated regime (stars), 
where we respectively have small and large depletions (also 
see Fig. [5] in Sect. 13.2(1 . The intermediate regime (dia- 
monds), where we have moderate depletions, is reasonable 
well described. 



The discrepancy in the thermally and sputtering dom- 
inated regimes can be fixed once it is realized that Udrift in 
Eq.|2Uis in fact the collecting speed of water molecules on 
the grain, v co iicct- In the thermally dominated or subsonic 
regime this collecting speed is set by the thermal motions 
of the gas. The grain is practically at rest with respect to 
the gas and only the thermal motions of the gas particles 
determine the rate at which ice is formed (this is expressed 
in Ea. 1121 see the case where sh 2 o <C 1 in this equation). 
In the sputtering dominated regime this collecting speed 
may effectively be set to zero since there is no net ice for- 
mation due to sputtering (i.e. no net collection of water 
molecules on the surface). In the intermediate regime the 
collection speed is simply the drift velocity of the grain. 
We thus rewrite Eq. into 



/(oo) = exp(— a(u 00 u BC t/u g as)njv(^o)/Qtrv) 



with 



^collect — 



max(« drift ,« th : 




H 2 0\ 



■ ^collect < «crit H2 

■ V collect > WcritH 2 



(31) 



(32) 



Here u C rit H IS ~9 km/s (see Sect. 13. 21) and we took « 
0.30 km/s. With these new equations the predictions by 
J&M are in better agreement with our model calculations. 
This is illustrated in the right panel of Fig. [21 where we 
have corrected for the effects of subsonic drift and sputter- 
ing. Notice that the thermally dominated regime and the 
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sputtering dominated regime are much better described 
now. 

In a part of the sputtering dominated regime, the pre- 
dictions of J&M still suffer from sputtering by He, CO 
and N2 at small values of A. In the intermediate regime 
they suffer from the onset of the thermal regime at the 
largest values of A. This leads to respectively somewhat 
larger and smaller values predicted by J&M compared to 
our models. The differences are very modest however. 

In order to calculate "drift m Eg. 1321 for the intermedi- 
ate regime, we may use 



"drift 



Mc 



(33) 



This equation is an approximation of Eq. 0] valid when 
"drift > "th- For a discussion on this equation and a deriva- 
tion o f it fr om first principles see for example IK wan fe Hi J 
l|l977|) and I Jura fe Morris! l|l985|) . 

In summary we conclude that in order to calculate 
/(oo) for a given (AGB) star, the original equation for 
/(oo) as derived by J&M will not suffice in many cases, 
unless thermal motions and sputtering are taken into ac- 
count. This can be done by using the modifications given 
in Eas.l31landl32l In the next section we will apply the re- 
sults of our model calculations, and the modified equation 
for /(oo) derived in this section, to individual stars. 

4. The model compared with observations : 
modeling individual stars 

In this section we will calculate the depletion of water 
ice from the circumstellar envelopes of different types of 
stars, based on Eas. 1271 and 1311 thru 1331 We will look at 
some of the details of the ice formation process around the 
modeled stars and discuss why ice does or does not grow. 
As examples we will consider an OH/IR star, a Red Super 
Giant (RSG) and a population II AGB star. 

The depletion, A, will be calculated in the following 
way. First, we determine "drift from Eg. 1331 and we check 
condition 1321 to obtain the correct "collect- Second, we cal- 
culate TTjv(r ) using Eq. |27| with w dust = "drift + " gas - 
Third, we use "collect and ttjv^o) to derive /(oo) in Eq. 
1311 Finally, we use Eq.[3n]to find A. All our models in this 
section will assume that do = 0.1 /im, p sl \ — 3300 kg m~ 3 , 
Qvv — 1, Qabs = 0.02 (here we follow J&M) , a — 1 and 
/d g = 0.01. 

For the condensation radius we take ro = 1600 AU in 
the case of the OH/IR and population II AGB star. This 
choice is based on the results of Fig. [3] and seems to be 
a r easonable value whe n compared with observations (see 
e.g. lSopka et al. ( 1985)). In the case of a RSG we will take 
ro = 16000 AU (see e.g. J&M and references therein). 



4.1. A typical OH/IR star 

A typical OH/IR star has a luminosity of L* = 10 4 Lq. 
It loses mass with rates between M = lO~ 5 M /yr and 



M = 10~ 4 M©/yr. The outflow speed of the gas is typi- 
cally i> gas = 15 km/s. Taking ro~1600 AU we find that the 
depletion of water vapour is given by A«5% and A«15% 
for M = 10~ 5 M Q /yr and M = 10~ 4 M o /yr respectively. 

We thus predict that a typical OH/IR star may show 
depletions of water vapour in its envelope up to at least 
15%. In the above calculations we find that the drift veloc- 
ity is in the order of "drift ~3 km/s or lower. This is near 
or in the intermediate regime where sputtering is unim- 
portant and the grain moves supersonically through the 
gas. The absence of sputtering allows the ice formation 
to occur. The predicted presence of ice is in qualitative 
agreement with the infrared spectra of many OH/IR stars, 
which indeed reveal the spectral bands of water ice with a 
variety of band strengths (see e.g. lSvlvester et ail l)l999(l '). 
Examples of typical OH/IR stars mode led in this section 
are AFGL 5379 (e.g.lOlivier et al.l200lh and OH 26.5+0.6 
(e.g. iMever et al.lll998|) . 

4.2. A typical Red Super Giant 

A RSG in some sense resembles the massive counterparts 
of AGB stars. Their MS mass is in the range 8 M Q <M < 
60 M Q ), and like AGB stars they are surrounded by a 
dusty envelope. It is interesting to see if ice can grow near 
these stars. 

A typical RSG star has a luminosity of L» = 
2xl0 5 LQ. It may lose mass with typical rates between 
M = 10~ 5 M o /yr and M = 10~ 4 M Q /yr. In extreme cases 
the mass-loss rate m ay near values of M = 10 _3 M Q /yr 
(see e.g. AFGL 4106 (|Molster et al.lll999il ). The outflow 
speed of the gas is typically w gas = 30 km/s. Taking 
ro~16000 AU we find that the depletion of water vapour 
is given by A = 0%, A«2% and A«8% for M = 
lO- 5 M /yr, M = 10- 4 M Q /yr and M = lO- 3 M /yr re- 
spectively. 

Based on the results above we predict that near RSGs 
hardly any water ice will grow. Typically we may expect 
A<2%. Only in extreme cases up to about 8% of the wa- 
ter vapout may be depleted into ice. The main reason for 
the inability of RSGs to form ice is that the drift velocity 
of the dust is generally too large. The above calculations 
show that for M<lO _4 M /yr the drift velocity will typi- 
cally be "drift~5 km/s or higher. This is in the sputtering 
dominated regime where sputtering heavily suppresses or 
completely blocks ice growth. The large drift velocities are 
mostly due to the high luminosities of RSGs. Only for ex- 
treme mass-loss rates, like M = lO _3 M /yr, efficient ice 
formation is possible since Udrift~l-5 km/s (intermediate 
regime) and hence sputtering is once more unimportant. 

The presence of water ice towards typical RSGs (e.g. 
a Ori and VY CMa) has to our current knowledge not 
been frequently reported in literature, altho ugh in the case 
of NM L Cyg water ice has been detected bv lMolster et alJ 
(2002). Post-RSGs with signs that they had high mass-loss 
rates when they where still in the RSG phase do indeed 
reveal the presence of ice (e.g. AFGL 4106, HD 179821 and 



C. Dijkstra et al.: Water ice growth around evolved stars 



13 



IRC+10420 iMolster et al.ll2002^ V This is in agreement 
with our findings in this section. 

4.3. A typical population II AGB star 

Finally we consider population II AGB stars which gener- 
ally have low metallicities. The ice formation in this case 
is limited by the amount of dust present in the outflow. 
The dust-to-gas ratio is at most equal to the metalicity, 
which can be as low as Z = 1 x 10~ 4 . Taking f^ g = 1CP 4 , 
r «1600 AU, w gas = 15 km/s, L* = 4xlO 3 L , and M = 
lO _4 M0/yr we find that the drift velocity of the dust is 
t> dri f t «0.5 km/s, which is close to the thermally dominated 
regime where ice formation is most efficient. However, the 
depletion of water vapour is only A = 0.1%, caused by the 
low dust-to-gas ratio. We thus find that the depletion in 
population II AGB stars will be low, even for high mass 
loss rate cases. 

5. Summary/Conclusions 

The results of our study can be summarised as follows. 

1. For a given grain size, ice formation becomes more effi- 
cient as the mass-loss rate of the central star increases. 
An increase in M increases the growth rate on the 
grains and reduces the effect of sputtering, promoting 
ice formation. 

2. In a circumstellar envelope, small dust grains will dom- 
inate the ice formation process if they are present in 
sufficiently large relative numbers. Compared to large 
dust particles, small particles have a better abillity to 
grow ice, since they are less sensitive for sputtering 
and, for a fixed M, have a larger total collecting area. 

3. For small particle sizes (ao = 0.01 /J.m), modest 
ice growth may already start at mass-loss rates of 
10 -6 M Q /yr and major ice formation only sets in when 
the mass-loss rates are typically of the order of several 
times 10~ 5 to 10 M^/yr, i.e. during the superwind 
phase of an AGB star. 

4. The ice formation process can be characterised in 
terms of thrift • For any given particle size, sputtering 
is an effective mechanism to heavily suppress or com- 
pletely block ice growth at thrift > 3 km/s. Between 
^drift^O.G km/s and 3 km/s ice formation becomes in- 
creasingly more efficient as Wdrift decreases, but no 
more than 20% of the water vapour will form ice. When 
^drift < 0.6 km/s, more than 20% of the vapour may 
condense. In this thermally dominated regime, it is M 
that controls the amount of ice that growths. 

5. Both crystalline and amorphous ice can be formed in 
circumstellar envelopes. In the outflow, initially crys- 
talline ice will form. When the dust temperature gets 
below ~65 K, the molecules are immobilised on the 
grain, preventing crystals from forming and resulting 
in amorphous ice. A substantial fraction of the grain 
may be composed of amorphous ice 66%). 



6. For an increasing mass-loss rate, ice formation starts 
closer to the central star, the final particle size will be 
reached earlier, and the transition point where amor- 
phous instead of crystalline ice will be formed moves 
inwards. In the meanwhile, the depth up to which UV 
photons from the interstellar medium can penetrate 
into the envelope, and make crystalline ice amorphous, 
moves outwards. An increase of the initial grain size 
generally has modest or no effect on the distance de- 
pendence of ice formation. Only the depth up to which 
UV photons from the interstellar medium can pene- 
trate into the envelope may change several orders of 
magnitude, increasing for larger particle sizes. 

7. For a star with a low optical depth envelope, the effect 
of the UV radiation is generally important throughout 
most of the envelope, and will most likely contribute 
to the formation of amorphous ice. However, the effect 
of the ISM can be cancelled in the inner regions, where 
the temperature exceeds 70 K. 

8. We have improve d upon an analytical equation, orig- 
inally derived bv l.lura fc Morrisl (^985), that predicts 
the depletion of water vapour from a circumstellar en- 
velope. 

9. Applying our model results and the modified equation 
of Jura & Morris to OH/IR stars, RSGs and popu- 
lation II AGB stars, we find that OH/IR stars reveal 
depletions up to at least 15%. Population II AGB show 
very little depletion due to the low dust to gas ratio 
of these stars. OH/IR stars and population II AGB 
stars mainly occupy the intermediate and the ther- 
mally dominated ice formation regimes. Water ice will 
not or hardly form in the circumstellar envelopes of 
RSGs due to the generally large drift velocities of the 
dust in their envelopes, caused by their high luminosi- 
ties. 
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